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© A method of synchronizing and channel estimation in a TDMA radio communication system. 

© The invention relates to a method of synchronizing a received signal frame in a TDMA radio system with a 
locally generated (TR) training sequence and of determining a channel estimate based on the comparison of the 
received signal frame and the locally generated training sequence. A first vector comprising M correlation values 
is calculated. In this vector the position m w of the centre of energy is determined. A second vector comprising N 
consecutive correlation values distributed around the centre of energy m w is chosen to form the channel 
estimate. The centre of energy in the second vector can be chosen as synchronizing position. 
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TECHNICAL FIELD . 

The present invention relates to a method of synchronizing and channel estimation in a TDMA radio 
communication system. * 

BACKGROUND OF THE INVENTION 

In TDMA radio communication systems (TDMA = Time Division Multiple Access) information is 
transmitted on a channel in the form of signal frames, that are transmitted by the transmitter during evenly 
distributed time intervals. In the spaces between these signat frames the transmitter is "silent". In order to 
synchronize the receiver to these signal frames each signal frame comprises a known synchronization word 
in predetermined positions within the signal frame. In for instance the European GSM system for mobile 
telephony this synchronization word is 26 bits long. When the receiver expects a new signal frame from the 
transmitter a training sequence, that is identical to the 16 central bits of the synchronization word, is 
generated by a training sequence generator in the receiver. The received signals are compared to the 
locally generated training sequence, and when the best possible correlation is obtained between this 
sequence and the received signals synchronization is considered to exist between the locally generated and 
the received signal. 

In addition to the synchronization itself the training sequence is also used for channel estimation; Since 
the radio channel often is subjected to multiple path propagation the receiver comprises some sort of 
equalizer to eliminate this phenomenon. The equalizer requires a time limited estimate of the impulse 
response of the channel. This impulse response can be obtained from the correlation signal. For this 
purpose not only a synchronization position but an interval or a window defining the multiple path 
propagation and indicating where the equalizer is to operate is required. 

A previously known method of finding the position of the window of the equalizer is to choose that 
interval of fixed length that contains most of the energy from the correlation. A drawback of this method is 
that disturbances in the received signal can give a correlation peak far away from the actual peak. For this 
reason the disturbance has a large instantaneous influence on the position of the window and creates a 
variance (uncertainty) in this position, which results in deteriorated receiver performance. 

SUMMARY OF THE INVENTION 

An object of the present invention is to provide a method of determining a more robust channel 
estimate in a TDMA radio communication system. 

This object is achieved by a method with the characterizing features of claim 1 . 

Another object of the invention is a method of determining the synchronization position in a signal frame 
of a TDMA radio communication system. 

This object is achieved by a method with the characterizing features of claim 10. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The invention, together with further objects and advantages thereof, may best be understood by making 
reference to the following description taken together with the accompanying drawings, in which: 

FIGURE 1 shows a correlation-time diagram illustrating a previously known method to determine the 

synchronization position and channel estimate in a TDMA radio communication system; 
FIGURE 2 shows a corresponding correlation-time diagram in a case where the radio channel is 

subjected to disturbance; : 
FIGURE 3 shows a correlation-time diagram illustrating the method of the invention in the disturbed 
case of Figure 2; 

FIGURE 4 shows a preferred embodiment of an apparatus for performing the method in accordance 
with the invention. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Figure 1 shows a correlation-time diagram, in which the sampling instances n run along the X-axis and 
the squared magnitude of the correlation between the locally generated training sequence and the received 
signal runs along the Y-axis. 
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The present invention has been described with reference to the European GSM-system. However, it is 
appreciated that the same principles also can be used in for instance the American mobile radio 
communication system in accordance with the standard IS-54. In fact the invention can be used in any 
TDMA-system that is based on synchronization of a received signal with a locally generated training 
5 sequence. - ■ 

It will be understood by those skilled in the art that various modifications arid changes may be made to 
the present invention without departure from the spirit and scope thereof, which is defined by the appended 
claims. 

io Claims 

1. A method in a TDMA radio communication system for choosing from a first vector comprising M 
correlation values between a synchronization sequence and M parts of a signal frame, which are 
partially overlapping and mutually displaced one sampling interval, a second vector of N correlation 
75 values for channel estimation, characterized by 

• - calculating the centre of energy w of said first vector in accordance with the formula 
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where c(k) designates correlation value k in said first vector; 

- rounding the calculated value w to the nearest integer for forming a preliminary window centre 
position m w in said first vector; and 

- choosing N consecutive correlation values distributed around said preliminary window centre 
30 position m w for forming said second vector. 

2. The method of claim 1 , characterized by 

- choosing 2L+1 partially overlapping vectors which are mutually time displaced one sampling 
interval, each comprising N = 2K + 1 consecutive correlation values distributed around said prelimi- 

35 nary window centre position m w ; and 

- choosing that of the 2L + 1 chosen vectors that has the largest energy content E„, that is that 
: vector that maximizes the expression 

K 



as said second vector. 

3- The method of claim 2, characterized by choosing L in dependence of the current time dispersion. 



4. The method of claim 3, characterized in choosing L proportional to the delay spread s defined in 
accordance with 

50. 



S = 
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The invention will be described below with reference to the European GSM-system. In this system a 
synchronization word comprises 26 bits. The 16 central bits in this word have good correlation properties 
when correlated with the entire synchronization word, that is a maximum correlation = 16 in the central 
position and a correlation of 0 in the remaining ten positions (C(k) = [0 0000 16 0000 6]). These 16 
central bits are generated as a training sequence in a training sequence generator in the receiver This 
training sequence is used for forming for instance 11 correlation values with the received signal frame. The 
squared magnitudes of these correlation values c(n) are shown in Figure 1. The finaJ synchronization 
position is chosen by comparing mutually displaced windows, each containing 5 correlation values, with 
respect to energy contents. The central position of that window that contains most energy is outputted as 
the synchronization position. Additionally the 5 correlation values c{n) within this window are outputted as a 
channel estimate to the equalizer. In Figure 1 the different windows have been indicated with reference 
designations 2-8, where these numbers indicate the central positions of the windows in the correlation 
vector that is formed by the 11 correlation values. In the example shown in Figure 1 the previously known 
method will therefore choose synchronization position m = 3, since the correlation peak lies at position 3 and 
the energy also is concentrated around this position. 

In practice the correlation peak is not as unambiguously defined as in Figure 1 . Rather the received 
signal isdisturbed by noise and by other transmitters, which results in disturbed channel estimates. Figure 
2 shows an unfavourable case, in which a disturbance peak has been introduced into correlation position 8. 
If this peak is sufficiently high the result of the previously known method can be that the total energy for 5 
correlation values happens to be larger around the disturbance than around the correlation peak. In this 
case the previously known method can therefore set m = 8 instead of the proper value 3. Furthermore, the 
channel equalizer will receive correlation values around the disturbance peak at m = 8, which means that the 
equalizer receives a very poor estimate, since the estimate is based on a disturbance peak. 

Figure 3 illustrates a preferred embodiment of the method in accordance with the invention. 

The method in accordance with the invention can be divided into several steps. In a first step the centre 
of energy w is calculated in accordance-with the formula: 

f>|c(*)| 2 



£ |c(Jc> |* 

where M is the number of correlation values. A suitable value for M has proven to be for instance the value 
11. The obtained value is rounded to a preliminary window position m w by rounding the obtained value w to 
the nearest integer. 

In a second step the energy of the correlation values c(n) that are contained in two windows around this 
preliminary central window position are calculated in accordance with the formula: 



k 

E n = £ |c(j + ^ +J2 )|2 ,3 = 0,1 

• J— * 

where 2K+.1 =N, that is the number of correlation values in each window, for instance 5. In Figure 3 this 
method implies that w will lie near 3, the preliminary window centre position will be rounded to 3, and two 
windows centered around positions 3 and 4 are compared with respect to energy contents. The coefficients 
c(n) of that window that has the largest energy content is outputted to the equalizer as a channel estimate. 

The final synchronization position m can be decided in several ways. One way is, as in the previously 
known method, to simply choose the centre position of the window. Another way is to determine the centre 
of energy in the chosen window in accordance with the formula: 
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where m w + n designates the central position in a chosen window. 

In the example described only two windows around the preliminary window centre position m* are 
w examined. However, it is possible to choose more windows than two. For instance 5 windows can be 
chosen around the preliminary window centre position m w . The number of windows can for instance be 
chosen in~ dependence of the current time dispersion. For instance the number of windows can be chosen 
proportionally to the delay spread s defined in accordance with the formula: 

J5 



s = 
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Furthermore it is appreciated that the number of correlation values in each window is not limited to 5. but 

25 that other values are also possible, 

An essential advantage of the present invention is that the determination of the centre of energy and the 
locking of the search window to values around this centre gives a more robust determination of both the 
synchronization position and the channel estimate. The reason for this is that a single disturbance peak in 
the correlation does not have such a strong influence on the calculation of the centre of energy. 

30 A preferred embodiment of an apparatus for performing the method in accordance with the invention 
will now be described with reference to Figure 4. 

The signal received by the receiver over an antenna is in a mixer converted to a base band signal and 
is in two A/D-converters divided into a complex input signal y(n) in a conventional way. This input signal is 
in a correlator CORR correlated with a training sequence from a training sequence generator TR local to 

35 said receiver. This training sequence generator outputs 16 values to correlator CORR. The correlation 
values c(n) of correlator CORR are stored in a buffer BUF1. In the example this buffer contains 11 memory 
positions for 1 T complex correlation values c(n). > 

In a unit 10 the square of the magnitude of each correlation value c(n) is formed. These squared values 
are outputted to a summation unit 20, which calculates a weighted average of the 11 input signals, where 

40 the weighting coefficients in integer steps run from 10 down to 0. Summation unit 20 can for instance be 
realized as an FIR-firter. In a similar way the sum of the squared magnitudes of the correlation coefficients 
c(n) are formed in a summation unit 30. This unit can also be realized as a FIR-firter, in which the filter 
coefficients are all set equal to 1. 

The output signal A from summation unit 20 is in a division unit 40 divided by the output signal B from 

45 summation unit 30 for forming the centre of energy w. In a rounding unit 50 w is rounded to the nearest 
integer m w . This value locks the reference position for a summation unit 60, for instance a FIR-firter, that 
performs the energy calculation for for instance 5 windows distributed around position m w and mutually 
displaced one sampling step, each window containing the squared magnitudes of each 5 correlation values 
c(n). These 5 sums, which represent the energy within each window, are stored in a buffer BUF2. In a unit 

5a MAX these 5 values are compared and the largest value is chosen, that is that window is chosen that had 
the largest energy content As output signal the centre position m in this window can for instance be 
outputted. In this embodiment this output signal also defines the synchronization position. Alternatively the 
centre of energy within the window can be calculated as indicated above for the final determination of the 
synchronization position. 

55 The window centre position m is also inputted to buffer BUF1, from which 5 correlation values c(n) 
centered around this position are chosen for transfer to a memory unit EST for storing the channel estimate. 
These 5 complex values are then transferred to the equalizer for channel equalization. 
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5. The method of claim 1, characterized by 

- choosing two partially overlapping vectors that are mutually time displaced one sampling intervaJ, 
each comprising N=2K+1 consecutive correlation values surrounding said preliminary window 
centre position m w ; and 

- choosing that of the two chosen vectors that has the largest energy content E„, that is that vector 
that maximizes the expression 



as said second vector. 



^\c(j+m w +n)\ 2 n = 0,1 



6. The method of any of claims 2-5, characterized by outputting the position of the centre position of said 
second vector in said first vector as final synchronization position m. 

7. The method of any of claims 2-5, characterized by 

- calculating the centre of energy x within said chosen second vector in accordance with" the 
formula - 

• ■ < K 



x = 

K 



J'-x - 

where m w +1 designates the centre position of said second vector in said first vector; 

- adding the value of x to the centre position m w + 1 of said second vector; and 

- outputting the position obtained in the previous step as final synchronization position m. 

8- The method of any of the proceeding claims, characterized by said first vector comprising 11 
correlation values, that is M = 1 1 . 

9. The method of any of the proceeding claims, characterized by said second vector comprising 5 
correlation vaJues, that is N = 5. 

10- A method in a TDMA radio communication system for determining from a vector comprising M 
correlation values between a synchronization sequence and M parts of a signal frame, which parts are 
partially overlapping and mutually time displaced one sampling interval, the synchronization position m 
of said fram in said vector, characterized by 

- calculating the centre of energy w in accordance with the formula 



M-l 



w = *=° 
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where c(k) designates correlation value k in said vector, said calculated centre of energy value w 
forming the sought synchronization position m. 
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